Canonizing curvature squared action in the presence of lapse function 

Abhik Kumar Sanyal 1 , Subhra Debnath 2 and Soumendranath Ruz 3 

August 31, 2011 



O 

£NJ ■ 1,2. Dept. of Physics, Jangipur College, Murshidabad, India - 742213 

3. Dept. of Physics, Satitara High School, Kandi, Murshidabad, India - 742170 

3 



< 



o 



oo 
o 



Abstract 



Lapse function appears as Lagrange multiplier in Einstein-Hilbert action and its variation leads to the (00) 
equation of Einstein, which corresponds to the Hamiltonian constraint equation. In higher order theory of 
gravity the situation is not that simple. Here, we take up curvature squared action being supplemented by 
an appropriate boundary term in the background of Robertson- Walker minisuperspace metric, and show how 
to identify the constraint equation and formulate the Hamiltonian without detailed constraint analysis. The 
5^ ' action is finally expressed in the canonical form, where, the lapse function appears as Lagrange multiplier. 

6J[), Canonical quantization yields Schrodinger like equation, with excellent features. 

PACS 04.50. +h 

> : 

1 Introduction 

OO ■ 

ly— ^ , Explaining the cosmic evolution, taking only geometric terms in the action, has turned out to be an important 
issue presently, since we do not have a scalar field at hand, as yet. Particularly, an action in the form, A = 
J[aR + j3R 2 + jR^ 1 ]^/— gd 4 x, apparently can challenge scalar field theories. The dominance of R^ 1 term at 
the very late stage of evolution leads to effective negative pressure, sufficient to explain the Snla data with 
t-H . an accelerating phase, the dominance of R term in the middle, keeps the nucleosynthesis and the growth of 
perturbation necessary for structure formation, unchanged from Friedmann model, while the dominance of R 2 
term in the early universe leads to inflation without invoking phase transition [1], [2]. Modification of the Einstein- 
Hilbert action by including curvature squared terms (R 2 , R^R^, C f _ l ^ p \C fJ - upX ), C^ yp \ being the Weyl tensor, is 
also important in many other respect. It leads to a renormalizable theory of gravity [3], even while interacting 
with matter [4] and is also asymptotically free [5] . Unitarity of higher derivative quantum theory of gravity has 
also been established [6]. Further, the Euclidean form of the Einstein-Hilbert action is not positive definite and 
therefore the functional integral corresponding to the ground state wave function of the universe diverges badly. 
A positive definite action that includes i? 2 -term, in the form [7], [8] 

1 



S = - l J <TX^ [ACf jkl +B(R- 4A) 2 ] , 

leads to a convergent integral for the ground state wave function, and it reduces to Einstein-Hilbert action in the 
weak field limit. Additionally, canonical quantization of the above action leads to a Schrodinger like equation, 
where an internal variable acts as the time parameter [7], [9], [10], [11]. A string inspired theory of gravity [12] 
and the 4-dimensional Brane world effective action [13] also contain such terms. In view of the above discussion, 
it turns out to be an important issue to include curvature squared term in the gravitational action and to study 
its quantum cosmological consequence, since it plays a dominating role only in the early universe. This requires 
canonical formulation of the theory. 

Hamiltonian formulation of a gravitational action containing curvature squared term is non trivial. Gibbons and 
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Hawking [14] have argued that the Einstcin-Hilbert gravitational action is incomplete, unless it is supplemented 
by an appropriate surface term that involves the integral over the boundary of space-time of the trace of the 
second fundamental form. To get a unique Hamiltonian of a theory, a unique boundary term is required which 
must be expressed as the integral over the boundary of space-time of some quantity. Without an appropriate 
boundary term the Hamiltonian and the corresponding quantum description of a theory usually turns out to be 
completely wrong [15]. The problem regarding canonical quantization of any other higher order theory of gravity 
is that, they are associated with different boundary terms [16] corresponding to different variational techniques, 
and apparently none of them is a total derivative term. Now, canonical formulation of any higher order theory of 
gravity also requires the introduction of an auxiliary variable. Sanyal [15] has shown that curvature squared action 
should be supplemented by a unique boundary term in the form, E = 4/3 j ( A R)K^fh d 3 x , where, K is the trace 
of the extrinsic curvature, h is the determinant of the metric of the three space and (3 is the coupling constant. 
Apparently, such a term is not a total derivative term. But then, E = o\ + U2 , where both can be individually 
expressed as total derivative terms. The first one, o\ — 4/3 j ( 3 R)K^/h d 3 x gets cancelled before the introduction 
of auxiliary variable and the other, viz., a 2 — 4/3 J ( 4 R — 3 R)K\fh d 3 x gets cancelled after the introduction of 
the auxiliary variable. This makes R 2 theory of gravity apparently as complete as Einstein-Hilbert action, when 
it is expressed as 



A = 13 J R 2 ^d 4 x + a 1 +a 2 . (1) 

Although, E has been found in view of Robertson- Walker minisuperspace metric, however it has also been tested 
successfully in a few other anisotropic minisuperspace metrics, viz., Kantowski-Sachs, axially symmetric Bianchi-I 
and Bianchi - III [15]. It is also important to note that the standard variational principle (metric formalism) 
of F(R) theory of gravity yields a surface term [17] in the form 2 J" E \fh F^Kd 3 x, which reduces to the above 
surface term E, for F(R) = R 2 and = 1. The surface term for F(R) theory of gravity also reduces to the 
Gibbons-Hawking term, when F(R) = lg ^ G . As already mentioned, E does not appear to be a total derivative 
term. What we have achieved in our earlier work [15], is to show that expressing E = cj\ + ct 2 and under suitable 
choice of auxiliary variable, the terms individually appear as total derivative terms. 

Still, there remains an unsolved important issue, which is presently our concern. The space-time metric ds 2 = 
g^vdx^dx" , under 3 + 1 decomposition can be expressed as, 

ds 2 = -(N 2 - N t N l ) dt 2 + 2N l dx l dt + h l3 dx l dx° , (2) 

where, N and Ni are the lapse function and shift vector respectively. In view of such decomposition, Einstein- 
Hilbert action when supplemented by Gibbons-Hawking surface term, 

^R-U^gJx+^fVKK**, (3) 

leads to a canonical action in terms of the basic variables hij and its canonical conjugate momenta irtj , in the 
form, 

A = J hij-K^ -H c - H ci d 3 xdt = J hij-K^ - NR. - N l Ui d 3 xdt, (4) 

where, the lapse function TV and the shift vector N l appear as Lagrange multipliers. Variation of the action (4) 
with respect to the shift vector gives the super momentum constraint, 

7U = 2DjTVi j = 0, (5) 

and variation with respect to the lapse function gives the super Hamiltonian constraint, 

H = (16ttG) G ijkmj ir kl - j^^h ( 3 R - 2A) = 0, (6) 

where, the metric on the superspace is expressed as, 

Gijki = -^—=(hikhji + huhjk — hijhki). (7) 

The canonical quantization is then straight forward, which gives the so called Wheeler-deWitt equation. Thus to 
claim that the action (1) is as complete as Einstein-Hilbert action (3), it should also be expressed in the canonical 
form (4). Such a construction has turned out to be impossible taking into account the whole superspace, and 
it has not so far been tried in minisuperspace model too. In the following section, we attempt this issue in the 
Robertson- Walker minisuperspace model, which accommodates lapse function only. In section 3, we extend our 
work to include Einstein-Hilbert action in addition to curvature squared term. 



2 Canonical formulation of (scalar) curvature squared action 

A general scale invariant action is expressed as, 

A = J [aC^ pX C^ px + 13 R 2 } V=gd 4 x. 

Since the Weyl tensor vanishes in the Robertson- Walker minisuperspace metric, 

dr 2 



ds 2 



-N{t) 2 dt 2 +a(ty 



1 — kr 2 



r 2 d8 2 +r 2 sin 2 ^ 2 



(8) 



hence the above action, in the presence of a cosmological constant A and supplemented by appropriate boundary 
term, reduces to 



g d x 



R 2 - 



2A 

J 



+ <Jl + 02, 



(9) 



The Ricci scalar, R^j^^ + ^z+N 2 -^- , under the choice hij = a 2 = z takes the form, 



R - w 



Z , T nfc ZN 

h A 2 

2z z 2zN 



(10) 



and the above action now reads, 
9z 2 



A = c/3 



36k 



Nz 



18Nzz 9N 2 z 2 36Nk 2 2NAz% 



dt + <7 1 +a 2 , (11) 



A 3 v/i \Ny/z~ N 2 ^J A 4 yi A 5 yi VI 

where, as already mentioned, 

E = a x + ct 2 = 4/3 J VhK 3 Rd 3 x + 4p J Vh K( 4 R - 3 R) d 3 x = 4/3 J Vh K 4 R d 3 x. 

and the constant c is the volume of the three space. Under integration by parts the first bracketed terms in the 
above action yield a counter term, that gets cancelled with cti and we are left with, 



(12) 



9z 2 l8Nzz 9N 2 z 2 18fci 2 36Afc 2 2AA^ 



dt + <7 2 ■ 



A 3 ^! A 4 VI A 5 VI Nz i yfi P 
where, we have chosen B — c/3 . At this stage we introduce the auxiliary variable, 
n _9A_ mR \ I Nz ' 

w dz [n 3 ^ a 4 vi ' 

and express the action in the canonical form as, 



A = B 



Qz Ai A 3 VI 2 18fcz 2 36Afc 2 2NAz% 



B BN® 36B 2 



-Q 



Nz? 



dt + (72 . 



(13) 



(14) 



(15) 



Now the first term is integrated by parts and the total derivative term gets cancelled with a 2 , and we are finally 
left with, 



• A 18Bkz 2 A 3 VI - /36Afc 2 2A Ar a 

-Qz - —zQ + —— 3 ^rQ + B — -= -Nz* 

A Nz? 36B \ y/z p 



dt. 



The canonical momenta are, 



A z 
p z = -Q - —Q + 36.Bfc— — 3-, p Q 

A Nz? 



■ Q 

-z, p N = -2 — . 



(16) 



(17) 
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The Q variation equation gives back the definition of Q given in (14). The z variation equation is, 



N 



—Q 

N 



iV2 

N 2 



Q - 36Bk 



Nzi 



Nz 
N 2 zl 



3z 2 



ANzi 



N 3 Q 2 
lIB^fl 



N 3B 



18Bk 2 — H — —NA^/z = 0, (18) 

Z 2 P 



and the N variation equation is 

zQ zQ z 2 ^N 2 

— — — + 18-Bk T + 

N N N 2 zi 



12B 



■Q 2 -B 



36k 2 



2A 3 



^rZ 2 = 0. 



(19) 



In view of the definition of momenta given in (17), it is clear that neither Q nor N is invertible, which signals 
the presence of a constraint in the theory. This is also apparent from that the fact that the Hessian determinant 
vanishes, i.e., H = | ^ ■ g ? ^ | = 0. Apparently, Q and N are having the same status in the action. Nevertheless, 
we know that Q is only an auxiliary variable and it has been introduced in the action keeping its canonical form 
intact. As a result Q variation equation gives back the definition (14) of Q and must not give any dynamics. 
Now, the third signal for the presence of constraint is that a particular variable is non-dynamical, i.e., none of 
the field equations contain second derivative of that variable. But then, we observe that second derivative of N 
appears in equation (18). This is definitely confusing. Finally, the fourth signal for the presence of a constraint is 
that one of the field equations must not be dynamical, i.e., it must not contain second derivative term. However, 
both the equations (18) and (19) contain second derivative terms. Thus, the action (16) contains a constraint, 
but the constraint equation is hidden. The reason is, unlike the case of Einstcin-Hilbert action, the action (16) 
contains first derivative of N. In fact, one can choose, a variable Q' = NQ to get rid of N term from the action 
(16). In the process, N acts as Lagrange multiplier and A^ term disappears from the z variation equation (18) 
while, z term disappears from the N variation equation (19). As a result, equation (19) being free from second 
derivative term, stands as the constraint equation. Instead, one can also remove z term from equation (19) in 
view of the definition of Q given in (14). Thus equation (19) takes the following form, viz., 



zQ 
' N 



NzQ 
N 2 



+ 18Bfc 



N 2 zi 



+ 



NWz 
36B 



Q 2 -B 



36k 2 

7^ 



2A a 
— z 2 



= 0. 



(20) 



This is the equation we were in search of, which does not contain second derivative term and hence is a constraint 
of the system under consideration. It can be easily verified that this is the Hamiltonian of the system in disguise, 



H r = N 



iQ NzQ z 2 



36k 2 

7^ 



2A s 
— z 2 



(21) 



which is constrained to vanish in view of equation (20). So, as in the case of Einstein-Hilbert action, here too the 
Hamiltonian can be obtained under the variation of the lapse function, but then one has to utilize the definition 
of auxiliary variable in addition. Now, the next question is, how to express the Hamiltonian in terms of the phase- 
space variables? This usually requires detailed constraint analysis [18], where, the constraint in the configuration 
space variable is equation (20) and that in phase space variable is Qpq — Np^ = 0, as is observed in view of the 
canonical momenta (17). However, we show that even without going into the details of constraint analysis, the 
Hamiltonian in terms of the phase space variables may be obtained in a straight forward manner. This is possible 
because N acts only as a Lagrange multiplier. The definitions of canonical momenta (17) yield, 

• N z 2 
PQPz = zQ + —zQ - 36Bk—^-. 

iv Nz 2 



So, 



-zQ- 



NzQ 
N 



18Bfc- 



Nzi 



-PqPz 



ISBk- 



Nz? 



-PqPz 



18Bfc 



where, we have replaced i by pq, instead of pn , since A^ acts as Lagrange multiplier and so the Hamiltonian 
must not contain p^ ■ Thus the Hamiltonian constraint equation in terms of the phase space variables is obtained 
as, 



H r = - 



PqPz 



Nz? 



36B 



36k 2 2A 2 
— = z 2 

sfz f3 



= 0. 



(22) 
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But still this is not the end of the story. To express the curvature squared action in the canonical form (4) as in 
the case of Einstein-Hilbert action, we need to express H c as H c = NH . But, first of all it is required to express 
the Hamiltonian in terms of the basic (configuration space) variables (instead of auxiliary variable) spanned by, 
{hij,TTij,Kij,ILij} which essentially are {z,p z , z,p z }, in the Robertson- Walker minisuperspace metric. For this 
purpose and to avoid confusion, the standard choice is z = x ( see [7], [9], [10], [11], [15]). Thus, our extended 
phase space is spanned by {z,p z ,x,p x }. But here, to express H c = NH, so that % remains free from N, we 
choose, 



z 

N' 



Hence, 



dA dAdx p x , . 
Q = ^ = didE^N^ and > V Q = -z = -Nx, 

So we need to replace Q by ^ and pq by — Nx in the above Hamiltonian. Thus, finally we are able to write, 
H c = N [x Pz + ^pl - —x- - 365^ + -j-z*) = NU = 0. (23) 

It is now straight forward to express the action (15) as, 

A = J (zp z + xp x - NH) dt d 3 x, (24) 

and we have achieved our goal of expressing the curvature squared action in the canonical form with respect to 
the basic variables. We can also anticipate a general form of the above canonical action. Remember, x — — 
2^2. = —2Kij, and so x — —1K^, where, is the extrinsic curvature tensor. Again, p x = — ^n^, being 
the momentum canonically conjugate to K^j . Hence we may write the above form of canonical action as, 

A = J (hijirv + kijU ij - NH) dt d 3 x, (25) 

where, 7r y is the momenta canonical to h, L j . This result that curvature squared action can be expressed in a 
general canonical form is of course new and exiting. Now, the canonical quantization of the Hamiltonian constraint 
equation is straight forward, which yields, 

in a* n 2 ( d 2 «a\ , nn , f x 2k\ T 2BAz T 

7i^ = -^(^ + ^J*- 18 ^U + -J* + ^- (26) 

Again under a further change of variable, the above equation takes the look of the Schrodinger equation, viz., 



9* fi z (\ d 2 n d 



da MB \ x dx 2 x 2 dx 



I2kx 24k 2 4A 2 

— — + — r _ Tir a3 

a 3 las opx 
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* = i? e * (27) 



where, n is the operator ordering index and a — z'^ = a 3 . Hence, the proper volume plays the role of internal 
time parameter. Note that the effective Hamiltonian 

h 2 ( 1 d 2 n d\ 

H e {x,a) = -—[-—^ + ^—)+V e (x,a), (28) 



(29) 



54B \x dx 2 x 2 dx, 
is hermitian, where the effective potential V e , given by. 



V e (x,a) = -B 



12kx 2Ak 2 4A 2 

— ~ + ~ ~ TJT a3 

a 3 xoe3 6px 



is a function of both the so called time variable a and x . The hermiticity of the effective Hamiltonian allows one 
to write the continuity equation for n = —1, as, 

^ + V.J = 0, 

oa 



G 



where, p 

ih 

1 



v]/*v]/ 



and J 



found only under the choice of the factor ordering index n 



(J x ,0, 0) are the probability density and the current density respectively, with, 
It is important to note that the continuity equation in the above standard form is 

-1 . Thus, factor ordering index has been fixed from 



physical argument. Further, there exists a continuous one parameter unitary group U(a) , which may be expressed 
by the Dyson series, Viz., 



U(a) = Texp 



H e (x, a)da 



where, T is the Dyson time ordering symbol. Finally, in the very early Universe when a is of the order of Planck's 
dimension, the term containing the cosmological constant remains subdominant in the effective potential and may 
be neglected. The extremization of the V e then yields, 



a = 



(t-to), 



(30) 



which has been obtained earlier [15]. This clearly depicts that power law inflation is an artefact of curvature 
squared action, and cosmological evolution must have started with a positive curvature parameter, k > 0. In the 
later epoch, A term starts playing a dominant role, but it is not possible to obtain a solution of the extremum of 
the potential in closed form, keeping the lambda term. 



3 Einstein-Hilbert action being modified by curvature squared term 

A general coordinate invariant fourth order action is expressed as, 
~R-2A 



I 



16ttG 



+ (3R J + aC„ upa C^ 



^gd i x. 



(31) 



This action can formally be cast to a unitary renomalizable quantum theory of gravity with positive energy states 
under Lee- Wick [5] prescription. In the Robertson- Walker minisuperspace metric the Weyl tensor vanishes and 
after being supplemented by appropriate boundary term, it is expressed as, 



R-2A 
IGttG 



+ m 2 



l~^gd X + G + o x + o- 2 , 



(32) 



where, a — j ^/hKd 3 x is the Gibbons-Hawking boundary term. This action leads to inflation without phase 
transition, followed by reheating [2]. Hence, if one is interested in canonical formulation of a general coordinate 
invariant fourth order action in Robertson- Walker minisuperspace metric, it is sufficient to start with action (32). 
It is important to note that action (32) may be cast into a positive definite one under appropriate choice of /3, 
and has a newtonian gravity long-distance limit as a classical theory. However, a more general quadratic action is 
expressed as [19], 



R-2A 
IGttG 



+ aR 2 + 1 R„„R» V + eR„ vpa R» vp ° + \UR 



l^gd x. 



(33) 



Instead of taking Kretschman scalar squared term, it is customary to express it as ex , where, 

x = 32^ / ( i W i ^ <7 - AR ^ V + r2 ) V=9d 4 x. 

X, also known as Gauss-Bonnet term is topologically invariant and so its functional derivative vanishes in four 
dimension. Further, OR is manifestly covariant total divergent term and is usually ignored. Hence, the above 
action reduces to, 



R-2A 
16ttG 



+ aR 2 + "/R^Br 



-gd A x, 



(34) 



which is also the action for the gravitational sector under BRST symmetry [20]. The above action may also be 
recast as, 



-gd x, 



(35) 
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just by choosing, j3 = a + A7. In the Robertson- Walker minisuperspace model, J(R^ l/ R 111 ' — \R 2 )^J —gd A x is 
again a total derivative term and thus is ignored. Hence, the above action (35) again reduces to the one given in 
(32), in the Robertson- Walker metric. It is interesting to note that the Weyl tensor squared term, 



R 



1 



:R 2 



may also be expressed as, 

C„ vpa C^° = (R^ pa R^»° - AR^ V + R 2 ) + 2(R^ - l -R 2 ), 

where, the first term is the Gauss-Bonnet term. Therefore, under the assumption of homogeneity and isotropy 
of space time, dubbed as the cosmological principle, both the actions (31) and (33) or (34) in disguise, reduce to 
(32), which we take up for the present purpose. In fact, instead of expressing (34) in the form given in (35), one 
could also express it as, 



A = 



D OA 1 



-gd x, 



which reduces to 

R- 2A 



16ttG 



+ PRuvW 



~gd x, 



in the Robertson- Walker minisuperspace model under consideration. This action is more interesting than (32), 
since in the gauge theoretic formulation of gravitation, R 2 only gives a scalar mode, while R^ V R^ U contains vector 
mode in the form of massless spin 2 gravitons. However, at present we pay our attention to action (32), which, 
as already mentioned, can be expressed as a positive definite action under appropriate choice of (3, and has a 
newtonian gravity long-distance limit as a classical theory. 

Following the same procedure as in section 2, one can express action (32) in the following canonical form, 



A = 



3i 2 



16ttG V 2Ny/z 



-6kNyG-2ANz% -Qz 



iVi iVVi 2 18Bkz 2 36BNk 2 



dt, (36) 



where, the auxiliary variable Q and the canonical momenta, pq and pn have the same expressions as found in 
equations (14) and (17) of the previous section, while p z is different, viz., 



3c / z 



Q 



yz 16ttG \Ny/z 
The TV variation equation now is, 
c / 3i 2 



AT 36Bkz 
—Q + — . 

Nzi 



WnG \2N 2 ^ 



6ky/z- 2Az5 



N N N 2 zi 



N 2 



12B 



■Q 



365fc 2 



0, 



(37) 



(38) 



which is supposed to be the Hamilton constraint equation. As before, removing z term, in view of the definition 
of the auxiliary variable (14), one can easily verify that this is again the Hamiltonian of the system in disguise, 



H, = N 



(- 



3i 2 



16ttG \2N 2 ^z 



6k^-2Az* - 



zQ NzQ 
~N ~ N 2 



\8Bk- 



N 2 zi 



+ 



365 



Q 2 - 



'MBk 2 

~7z~ 



(39) 



which is constrained to vanish. Now, using the expression, 



PQPz 



+ 



3cz 2 



N 



36£fci 2 

~a^I~ : 



16wGN^ N ^ 

the Hamiltonian constraint equation in terms of the phase space variables is obtained as, 



16ttG I 2Ny/l 



r- A 18-Bfc 2 

SNkyfz + 2NKz* - pqp. 



Nz? 



pq 



3673 W 



mBNk 2 



0. 



(40) 



(41) 
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Finally, to express H c = NW. , let us choose as before, 



z 

X= N 



Q = 



dA dA dx p x 



dz 



dxdz = N' and > PQ = -* = - N *> 



So we need to replace Q by ^ and pq by — Nx in the above Hamiltonian. Thus we get, 



H r = N 



xp z + 



Vz_ 
36B 



Px - 



WBk 



x -36B— + 



3x 2 



Jz 16ttG \2y/z 



- 6ky/z + 2Az 



NH = 0. 



(42) 



It is important to notice that momentum is not associated with the Einstcin-Hilbert sector of the action. The 
action (36) can now be expressed in the canonical form with respect to the basic variables as, 



= /(ip, 



xp x - NH) dt d 3 x = 



The corresponding quantum version is 



J~z~ dz 



36Bx 



dx 2 



d 



-—)*-18Bk[ — + — 

x ox 



2k 



NH) dt d 3 x 



(43) 



16ttG 



3x 
~2z~ 



k . z 
6- +2A- 

x x 



(44) 



Again under a further change of variable, the above equation takes the look of the Schrodinger equation, viz. 



if l 



da 



54B 



1 d 2 

x dx 2 



d_ 

dx 



16ttG 



4fc 4Aa3 
x 3x 



- \2Bk 



x 

T 

a 3 



2k 



XCt3 



* = ff e *,(45) 



where, a = z? again plays the role of internal time parameter and the effective potential V e , given by, 



V e 




a 3 XOL3 



(46) 



is a function of both the so called time variable a and x, as before. We observe that Einstein- Hilbert sector 
appears as effective potential in the quantum version of the theory. We can try to explore the reason for this. In 
the case of E-H action alone, the momentum is p 2 oc i. But here, under the introduction of auxiliary variable, 
one gets, pq = —z, while a complicated expression for p z , gives Q. With the choice of basic variable, this pq has 
been transformed to a generalized co-ordinate x , which finally appears in the effective potential. 
The effective Hamiltonian (H e ) is again hermitian and so, one to write the continuity equation for n = — 1 as 
before, and the probability interpretation follows. 



4 Summary 

Summarily, in an earlier work [15], we have presented a unique boundary term to supplement curvature squared 
action like the Gibbons-Hawking term, which supplements Einstein-Hilbert action. This boundary term is at par 
with the one obtained under the variation (metric formalism) of F(R) theory of gravity, choosing F(R) = R 2 . 
Such a boundary term usually can not be expressed as a total time derivative term. However, we have shown that, 
the boundary term X = o\ + a 2 , where, both are total derivative terms, under appropriate choice of auxiliary 
variable. In view of such a boundary term, here, we have been able to express the curvature squared action 
in the canonical form, in the presence of lapse function, taking Robertson- Walker minisuperspace model under 
consideration. The canonical form of the action is in close resemblance with the one for Einstcin-Hilbert action. 
The phase space of the curvature squared action, being a higher derivative theory of gravity, is larger and is 
spanned by hij ■, K i3 ;, it 13 W 3 ' . Canonical quantization then leads to a Schddinger like equation, where, a = a 3 , 
the proper volume, acts as the internal time parameter. The effective Hamiltonian is self adjoint and so the 
probabilistic interpretation is straight forward. To express the continuity equation in the standard form, one has 
to fix up the operator ordering parameter n = — 1. Since the effective Hamiltonian has time dependance through 
the effective potential, so the unitary transformation is given by the Dyson series. Finally, the extremization of 
the effective potential gives coasting solution for the curvature parameter k = 1 , which implies that inflation is 
an artefact of curvature squared action for a initially closed Universe. Although, the whole process of canonizing 
has been carried out in a minisuperspace model, however, it might give insight to extend the work in the whole 
supcrspace. 
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